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Abstract 

Background: The variation in tine rate of intraoperative respiratory events (IRE) over tinne under anesthesia and the 
influence of anesthesia-related factors have not yet been described. The objectives of this study were to describe 
the risk over time and the risk factors for IRE in children at a tertiary care hospital in southern Thailand. 

Methods: The surveillance anesthetic database and chart review of IRE of 14,153 children who received surgery at 
Songklanagarind Hospital during January 2005 to Decennber 201 1 were used to obtain demographic, surgical and 
anesthesia-related data. Incidence density of IRE per person-time was determined by a Poisson modelling. Risk of 
IRE over time was displayed using Kaplan Meier survival and Nelson-Aalen curves. Multivariate Cox regression was 
employed to identify independent predictors for IRE. Adjusted hazard ratios (HR) and their 95% confidence intervals 
(CI) were obtained from the final Cox model. 

Results: Overall, IRE occurred in 315 out of 14,153 children. The number (%) of desaturation, wheezing or 
bronchospasm, laryngospasm, reintubation and upper airway obstruction were 235 (54%), 101 (23%), 75 (17%), 21 
(5%) and 4 (1%) out of 315 IRE, respectively. The incidence density per 100,000 person-minutes of IRE at the 
induction period (61.3) was higher than that in the maintenance (13.7) and emergence periods (16.5) (p< 0.001). 
The risk of desaturation, wheezing and laryngospasm was highest during the first 1 5, 20 and 30 minutes of anesthesia, 
respectively. After adjusting for age, history of respiratory disease and American Society of Anesthesiologist (ASA) 
classification, anesthesia-related risk factors for laryngospasm were assisted ventilation via facemask (HR: 18.1, 95% CI: 
6.4-51.4) or laryngeal mask airway (HR: 12.5, 95% CI: 4.6-33.9) compared to controlled ventilation via endotracheal tube 
(p < 0.001), and desflurane (HR: 1 1.0, 95% CI: 5.1-23.9) compared to sevoflurane anesthesia (p < 0.001). 

Conclusions: IRE risk was highest in the induction and early maintenance period. Assisted ventilation via facemask or 
LMA and desflurane anesthesia were anesthesia-related risk factors for laryngospasm. Therefore, anesthesiologists 
should pay more attention during the induction and early maintenance period especially when certain airway 
devices incorporated with assisted ventilation or desflurane are used. 
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Background 

Despite the low rate of mortality associated with cardiac 
arrest in anesthetized children over the past few decades 
[1], as many as 27% of perioperative cardiac arrests were 
caused by respiratory events in the Pediatric Perioperative 
Cardiac Arrest registry [2]. Tait et al. [3] reported an 
incidence of 1.7% for respiratory events during the in- 
traoperative period in a retrospective study in children 
with upper respiratory tract infection (URI). A Thai 
Anesthesia Incidents (THAI) prospective study also 
reported a similarly low incidence (1.5%) of perioperative 
respiratory events (PRE) in 25,098 Thai children [4]. In 
contrast, four recent prospective studies have reported 
incidences of 15% for PRE in general child populations 
or as high as 28% in high risk children such as those 
with URI or obesity [5-8]. However, the separate incidence 
of intraoperative respiratory events (IRE) was mentioned 
in only one of those studies [6]. Known risk factors for 
PRE include history of upper respiratory tract infection, 
age less than 6 years, obesity, airway surgery, and absence 
of a specialized pediatric anesthesiologist during the op- 
eration [5-8]. However, descriptive profiles concerning 
changes in the rate of IRE over time in children as well 
as incidence density of IRE have never been described. 
Parents often ask whether the anesthetic risk increases 
with time under anesthesia. The period of anesthesia 
that anesthesiologists should be most concerned about 
is still debatable. Moreover, the predictive ability of some 
characteristics of anesthesia, such as anesthetic technique, 
related to the change of rate of IRE over time, has also 
never been investigated. Therefore, the objectives of this 
study were to describe the risk profile of IRE over time 
and identify the risk factors for IRE related to patient, sur- 
gical and anesthesia profiles. 

Methods 

A historical cohort study was conducted at Songklanagarind 
Hospital, an 853-bed tertiary care hospital in southern 
Thailand, after approval by the Ethics Committee, Prince 
of Songkla University [EC 552100813]. Written informed 
consent to participate in the study was waived. The sur- 
veillance anesthetic database under quality assurance by 
nurse anesthetists was used to identify IRE in children 
aged less than 15 years undergoing surgical procedures as 
well as digital imaging or other interventions under general 
anesthesia (GA) with or without regional anesthesia or 
peripheral nerve block between the period of January 
2005 and December 2011. All identified IRE were con- 
firmed by chart review of either paper-based or electronic 
anesthetic records which were available since 2005. The 
anesthetic records are divided into 2 sections; intraopera- 
tive section and post-anesthetic care unit (PACU) section 
where vital signs are routinely recorded every 5 minutes 
thoughout the anesthetic procedure. The timing of any 



respiratory events or rescue medication is recorded in 
the event section in the anesthetic record. 

Participants 

All procedures including regional anesthesia or peripheral 
nerve block were done by a certified anesthesiologist with 
at least one year s experience. Patients were excluded if they 
were classified as American Society of Anesthesiologists 
(ASA) physical status 4 or 5, had a preoperative arterial 
oxygen saturation at room air < 95%, were already endo- 
tracheally intubated and/or mechanically ventilated prior 
to surgery or had congenital heart disease with cyanosis. 
These criteria meant that existing hypoxemic children 
or children who previously had severe respiratory or 
circulatory problems were excluded. 

Anesthesia practice 

In routine practice in the hospital during the period of 
study, the choices of anesthetic technique, airway devices 
and anesthetic agents were at the discretion of the at- 
tending anesthesiologists. Attending anesthesiologists 
supervised all anesthesia residents and nurse anesthetists 
for all anesthetized children from the start of anesthesia 
until discharge from the PACU. There was a total of 28 
attending anesthesiologists throughout the 7 year study 
period. The average number of attending anesthesiologists 
each year was 17; 2 were certified pediatric anesthesiol- 
ogists, 8 had experience in pediatric anesthesia > 5 years 
and 7 had experience < 5 years. The induction technique 
for GA was either inhalation using sevoflurane or intraven- 
ous injection. The maintenance of anesthesia was either by 
spontaneous breathing with 100% oxygen using intravenous 
anesthesia, or by facemask or laryngeal mask airway (LMA) 
using inhalation anesthetic agent with N2O-O2 or Air-02, 
or by endotracheal tube (ETT) using neuromuscular block- 
ing agent (NMBA). 

Standard operating procedures 

The mode of ventilation was either by manually assisted 
ventilation or controlled ventilation by anesthetic machine 
ventilator. For manually assisted ventilation via facemask 
or LMA, tidal volume was adjusted at 7-10 ml/kg by 
anesthesia personnel and synchronised with children's 
spontaneous breathing. Assisted ventilation by anesthetic 
machine ventilator was not applied if children could 
breathe spontaneously. For ETT under controlled mechan- 
ical ventilation, the anesthetic machine ventilator was set at 
the same tidal volume. Positive end-expiratory pressure was 
not routinely applied during controlled ventilation in chil- 
dren except in prolonged (more than 3 hours) anesthesia. 

One hundred percent oxygen with high flow (5 liters per 
minute) was applied at least one minute before intubation 
or LMA insertion and immediately after extubation. Basic 
intraoperative monitoring including pulse oximetry and 
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end tidal carbon dioxide were used during the period of 
anesthesia. After the operation, children were sent either to 
the PACU or to the pediatric intensive care unit (PICU) de- 
pending on the decision of the surgeon or anesthesiologist 
One hundred percent oxygen was not usually administered 
while relocating children to the PACU unless persistent de- 
saturation (arterial oxygen saturation < 95% with 100% oxy- 
gen more than one minute) occurred before transportation. 

Outcome of interest and variables 

Laryngospasm, wheezing or bronchospasm, hypoxemia 
or desaturation (oxygen saturation < 95% for more than 
10 seconds [9]), upper airway obstruction (UAO; soft 
tissue obstruction or secretion obstruction after the airway 
device was removed) and re-intubation were chosen 
as significant events constituting IRE based on their 
seriousness which can lead to severe hypoxemia and 
cardiac arrest in children. The categories of IRE were 
not mutually exclusive. Thus, laryngospasm, bronchospasm 
and UAO might have desaturation as a result. Time to first 
IRE was the outcome of interest, assessed by the researcher 
by performing chart review of all IRE anesthetic records 
during the data management period. Cause of desaturation 
was also assessed. 

Variables collected were categorised into patient-related 
risk factors, surgical-related risk factors and anesthesia- 
related risk factors. Patient-related risk factors were age, 
sex, weight (kg), height (cm), body mass index (kg/m^), 
recent upper respiratory tract infection (URI), pulmonary 
disease such as asthma or pneumonia, allergic rhinitis, 
food or drug allergy, non-cyanotic heart disease, anemia 
(hematocrit < 39% in male and < 36% in female). Surgical- 
related risk factors were type of surgery (elective/ 
emergency), type of patient (inpatient/outpatient), position 
(supine, lateral, lithotomy, prone) and site of operation 
(superficial; chest wall, abdominal wall and minor surgery, 
bronchoscopy or laryngoscopy, ear-nose-throat, eye, 
intrathoracic, intraabdominal, spinal or extremity, re- 
mote; radiology services and cardiac catheterization). 
Anesthesia-related risk factors were ASA classification, 
choice of anesthesia (GA, GA with regional anesthesia, 
GA with peripheral nerve block), premedication agent 
(none, chloral hydrate, diazepam, others), induction 
anesthetic agent (propofol, thiopental, sevoflurane, others), 
intubating agent (non-depolarizing muscle relaxant, 
succinylcholine, sevoflurane), narcotic (fentanyl, morphine, 
others, none), airway device with mode of ventilation 
(ETT-controlled ventilation, jet ventilation-controlled 
ventflation, facemask-assisted ventilation, LMA-assisted 
ventilation), inhalation agent (sevoflurane, isoflurane, 
halothane, desflurane), gas mixed with oxygen (Air-02, 
N2O-O2, 100% O2), use of a NMBA and experience of key 
attending anesthesiologist (< 5 years, > 5 years, certified 
pediatric anesthesiologist). 



Missing data from the anesthetic database were completed 
by hospital information system or the electronic anesthetic 
record by the researcher. Height was missing fi:om 1.4% of 
records and was replaced with the mean height of the same 
age and sex. 

Definitions of periods of anesthesia 

The period of anesthesia of the 14,153 records and the 
timings of first IRE was divided into 3 periods of anesthesia; 
induction, maintenance and emergence periods. The in- 
duction period was defined as the time between starting 
of anesthesia and the starting of incision. The maintenance 
period was defined as the time between starting of incision 
and finishing the operation. The emergence period was 
defined as the time between finishing the operation 
and the airway device being removed. The 3 periods of 
anesthesia were defined using variable cut points of 
time based on the total duration of anesthesia adjustable 
by airway device used. For the total duration of anesthesia 
more than 30 minutes, the induction and emergence 
periods were 15 minutes after starting anesthesia and 
15 minutes before ending of anesthesia, respectively, 
because most chfldren received ETT intubation or LMA 
insertion. For children in which the total duration of 
anesthesia was less than 30 minutes, the induction and 
emergence periods were 5 minutes after starting anesthesia 
and 10 minutes before ending anesthesia, respectively, 
because most children received facemask ventilation 
with spontaneous breathing. The maintenance period 
was calculated based on subtraction from the total period 
of anesthesia of the combined duration of induction and 
emergence periods. 

Statistical analysis 

Descriptive statistics were used to summarize character- 
istics of the patients. The incidence density of any IRE 
per person-time (person-minutes) based on period of 
anesthesia and duration of each period was calculated 
by fitting a Poisson regression model. Significant differ- 
ences of incidence density of any IRE among the periods 
of anesthesia were determined if the p-value from the 
likeUhood ratio test was less than 0.05. 

Time to first occurrence of IRE was calculated from 
start of GA until IRE or, if no IRE occurred, patients were 
considered as censored at the end of GA. Risk of occur- 
rence of IRE over time was displayed using Kaplan Meier 
and Nelson- Aalen curves [10]. Log rank tests were initially 
employed to assess the differences in time-to-event across 
the subgroups of potential predictor variables. Cox re- 
gression models were fit to identify independent predic- 
tors of IRE based on the timings of first events of any IRE, 
laryngospasm, desaturation, wheezing or bronchospasm, 
or UAO in combination with reintubation. After run- 
ning the multivariate Cox regression model, variables 
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which violated the proportional hazard assumption were 
considered to have time-varying effects and explored 
further to determine if their effects differed among the 3 
periods (induction, maintenance and emergence). This was 
done by creating 3 partial proxy variables for each variable 
violating the proportional hazard assumption and fitting 
them in place of the original variable. Adjusted hazard 
ratios and their 95% confidence intervals (CI) were ob- 
tained from the final Cox model. In multivariate mod- 
eling techniques, an association with the outcome was 
considered statistically significant if the likelihood ratio 
p-value was less than 0.05. 

Sample size considerations 

Adequacy of sample size for identifying the risk factors for 
IRE was based on a preliminary examination of a subset of 
the children. The IRE-free Kaplan-Meier profile of 2,109 
children undergoing GA in 2010, among whom 44 events 
were recorded, revealed an event-free probability of ap- 
proximately 0.975 corresponding to a GA duration of 
2.5 hours. With a total sample of approximately 14,000 
children and a prevalence of exposure to a risk factor 
between 5% and 50%, differences in 1-hour risk of IRE 
as small as 2% to 1% (depending on prevalence) should 
be detectable as significant at an alpha of 0.05 and a 
power of 80%. We therefore believed that a sample size 
of 14,000 would be adequate. 

Results 

A total of 16,316 children undergoing GA from 2005 to 
2011 were identified under the period of data management 
between May 2012 and December 2012. We excluded 
1,552 cyanotic heart disease and/or required endotra- 
cheally intubated prior to surgery, 507 ASA classifica- 
tion > 4 and 104 preoperative arterial oxygen saturation at 
room air < 95%. IRE occurred in 315 (2.23%) of the 14,153 
study children, the numbers (%) of desaturation, wheezing 
or bronchospasm, laryngospasm, reintubation and UAO 
events were 235 (54%), 101(23%), 75 (17%), 21 (5%) and 4 
(1%), respectively. The causes of intraoperative desaturation 
(235 children) were light anesthesia (children were not 
adequately anesthetized) with endotracheal tube placed 
(31%), laryngospasm (26%), bronchospasm (22%), prolonged 
apnea during application of laryngoscopy (12%), accidental 
extubation (4%), endobronchial intubation (1.4%) and UAO 
(1.4%). The median duration of anesthesia was 90 minutes 
(interquatile range (60,140 minutes)). One hundred and 
nineteen children developed more than 1 IRE. Table 1 
shows the number and incidence density of IRE per 
100,000 person-minutes among each period of anesthesia. 
The incidence density of each type of IRE except UAO in 
combination with reintubation was highest in the induction 
period of anesthesia (p < 0.001). 



Table 1 Number and Incidence density of intraoperative 
respiratory events (IRE) per 100,000 person-minutes 
among 14,153 children 



Period of anesthesia 


Induction 


Maintenance 


Emergence 


p-value 


(person-minutes) 


(200,584) 


(1,156,812) 


(206,453) 




IRE 


Number (incidence per 100,000 






person-minutes) 






Desaturation 


90 (44.9) 


124 (10.7) 


21 (10.2) 


< 0.001 


Wheezing 


48 (23.9) 


41 (3.5) 


12 (5.8) 


< 0.001 


Laryngospasm 


27 (13.5) 


38 (3.3) 


10 (4.8) 


< 0.001 


UAO & reintubation 


1 (0.5) 


15 (1.3) 


7 (3.4) 


0.056 


Any 


123 (61.3) 


158 (13.7) 


34 (16.5) 


< 0.001 



p-value from likelihood ratio test. 
UAO: Upper airway obstruction. 



Characteristic profiles of intraoperative desaturation, 
wheezing, laryngospasm and UAO in combination with 
reintubation 

Figure 1 shows the IRE-free Kaplan-Meier curve. The 
probabiUty of being IRE free after 500 minutes was 0.92. 
IRE was found in all periods of anesthesia including the 
induction (38%), maintenance (43%) and emergence 
periods (19%). However, these percentages do not reflect 
the actual risk of IRE during these periods since the 
number of subjects at risk in the later periods declines 
over time. More appropriately. Figure 2 shows the Nelson 
Aalen plot of cumulative hazard of IRE describing the 
change in rate of different forms of IRE over time by the 
slope of the cumulative curve. 

The curve for desaturation was steepest during the 
first 15 minutes, while upward convexity was evident in 
the first 100 minutes indicating lowering of risks over 
this period. After 100 minutes, the convexity gradually 
disappeared suggesting that the risk had stabilized. 

Other IRE included wheezing, laryngospasm and UAO 
and reintubation. Their cumulative hazard profile had 
different shapes. The profile for wheezing was steepest 
during the first 20 minutes and stabilized after 1 hour. After 
300 minutes, the curve had somewhat of an upward con- 
cave shape, suggesting increasing risk after around 5 hours 
of GA. The curve for laryngospasm was steepest during 
the first 30 minutes and then stabilized with no laryngos- 
pasm occurring after 400 minutes. By contrast, the curve 
for UAO in combination with reintubation exhibited an 
upward concave shape, indicating a steadily increasing risk 
over the first 180 minutes; however, there was no subse- 
quent event after that time. 

Univariate and Cox regression analysis of IRE 

The univariate Kaplan Meier analysis for overall IRE 
included patient-related (8 variables), surgery-related 
(4 variables) and anesthesia-related (11 variables) risk 
factors (see Additional file 1). The overall IRE probability 
(95% CI) of IRE within 60 minutes of anesthesia was 0.017 
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Figure 1 Kaplan Meier survival curve of the first intraoperative 
respiratory event. 



(0.015, 0.019). Fourteen potential predictors which had log 
rank p values < 0.2 were included in the initial multivariate 
Cox regression analysis for IRE. After model refinement, 
8 out of the 14 potential predictors remained in the final 
model (see Additional file 1). Age, history of respiratory 
disease and ASA classification violated the assumption 
of proportional hazards and each was therefore used to 
create 3 partial proxy variables to allow their effects to differ 
among induction, maintenance and emergence periods and 
the dataset reorganized into multirecord format. Table 2 
shows the results of multivariate Cox regression of IRE 



n at risk 0 

Desat 14153 
Wheezing 14153 
Laryngo 14153 
UAO & Rein 14153 




Figure 2 Nelson-Aalen plot of the first event of different forms 
of intraoperative respiratory event. Desat; Desaturation, Laryngo; 
Laryngospasm, UAO & rein; Upper airway obstruction in 
combination witli reintubation. 



incorporating these time-varying effects. The significant 
risk factors were lateral position, airway surgery, airway 
device with mode of ventilation, use of succinylcholine 
as an intubating agent and desflurane anesthesia, whereas 
age, history of respiratory disease and ASA classification 
were allowed to have time-varying effects. It should be 
noted that children whose airway received controlled 
ventilation via ETT had 2.0- to 2.4-fold lower risk than 
other groups (p < 0.001). During induction period, an 
elevated risk was associated with age < 1 years (p = 0.004), 
having pulmonary disease (p < 0.001) and ASA classification 
3 (p < 0.001). In the maintenance period, risk was elevated 
in patients age < 6 years (p < 0.001) and having URI 
(p = 0.008), while during the emergence period, a higher 
risk was associated with age < 1 years (p = 0.048). 

Tables 3, 4, 5, 6 show the risk factors for different forms 
of IRE. The risk factors for desaturation were similar to 
those for any IRE. The risk factors for wheezing were 
having URI or pulmonary disease (p < 0.001) and ASA 
classification 3 (p = 0.014). The risk factors for laryngospasm 
were having pulmonary disease (p = 0.02), assisted venti- 
lation via facemask or LMA compared to controlled 
ventilation by ETT (p < 0.001) and desflurane anesthesia 
(p < 0.001), whereas age and intubating agent were the 
factors with time-varying effects. The probability of lar- 
yngospasm during the induction period was increased 
by the use of succinylcholine (p < 0.001), in the mainten- 
ance period by 1-6 years of age (p = 0.004) and at the 
emergence period by < 6 years of age (p = 0.002). The only 
risk factor for UAO in combination with reintubation was 
airway surgery, whereas age was an effect-varying factor 
with probability of UAO in combination with reintuba- 
tion in the maintenance and emergence periods increased 
by < 1 years of age (p = 0.01 and p = 0.005, respectively). 

Discussion 

The overall proportion of IRE in children undergo- 
ing GA between 2005 and 2011 in Songklanagarind 
Hospital was 2.23%, respectively. A longer anesthesia time 
(from 1 to 8.7 hr) correlated with a higher IRE probability 
(from 1.7% to 8%). The low incidence of IRE compared to 
that of other studies could be explained by the different 
criteria used to define desaturation in previous studies 
[4,8,11]. Although the proportion of IRE was highest at 
the maintenance period of anesthesia (43%), the incidence 
density per 100,000 person-minutes of any IRE was high- 
est (61.3) at the induction period compared to the main- 
tenance (13.7) and emergence periods (16.5). 

The patient and surgery-related risk factors for IRE, 
namely age < 6 years, having respiratory disease and 
undergoing airway surgery were quite similar to those 
reported in other studies [4-6,8]. However, history of aller- 
gic rhinitis, history of food or drug allergy and experience 
of key anesthesiologist were not found to be associated 
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Table 2 Multivariate Cox regression of IRE based on period of anesthesia (n = 42,055) 



Variable 


All periods (315) 


Induction (123) 


Maintenance (158) 


Emergence (34) 




Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Age (ref= > 6) 




1' 0.004 


l'< 0.001 


1' 0.048 


1-6 




1 .5 (0.94, 2.3)' 


2.3 (1.6, 3.4)^ 


2.0 (1.0, 4.1)' 


< 1 




2.5 (1.5, 4.2)^ 


3.0 (1.8, 4.8)^ 


3.1 (1.2, 8.2)^ 


Respiratory disease (ref = no) 




l'< 0.001 


1' 0.008 


1' 0.2 


URI 




1 .4 (0.67, 2.9)' 


2.2 (1.3, 3.7)^ 


2.1 (0.72, 6.3)' 


Pulmonary disease 




3.5 (2.3, 5.5)^ 


1.6 (0.9, 2.8)' 


2.1 (0.81, 5.7)' 


Site of operation (ref = superficial) 


l'< 0.001 








Broncho/laryngoscopy 


3.6 (2.0, 6.6)^ 








Ear-nose-throat 


1 .6 (0.86, 2.8)' 








Eye 


1 .0 (0.55, 1 .8)' 








Thoracic 


1.1 (0.59, 2.2)' 








Abdomen 


1 .4 (0.83, 2.4)' 








Spine & extremity 


1 .3 (0.74, 2.3)' 








Remote* 


1 .3 (0.48, 2.8)' 








Position (ref = supine) 


1' 0.009 








Lateral 


2.0 (1.3, 3.1)^ 








Lithotomy 


0.50 (0.23, 1.1)' 








Prone 


1.1 (0.51, 2.2)' ^ 








ASA classification (ref= 1) 




l'< 0.001 


1' 0.83 


1' 0.75 


2 




1.7 (0.96, 3.1)' 


1.0 (0.66,1.5)' 


1 .3 (0.64, 2.7)' 


3 




4.0 (2.1, 7.8)^ 


1.1 (0.66, 2.0)' 


1 .2 (0.38, 3.8)' 


Airway device (ref = E^-C) 


l'< 0.001 








Jet ventilation-C 


2.4 (1.1, 5.1)^ 








Facemask-A 


2.4(1.4, 4.0)^ 








LMA-A 


2.0 (1.2, 3.3)^ 








Intubating agent (ref=NDMR) 


l'< 0.001 








Succinylcholine 


1.7 (1.2, 2.2)^ 








Sevoflurane 


1 .6 (0.9, 2.7)'^ 








Inhalation agent (ref = sevoflurane) 


1' 0.0013 








Isoflurane 


0.77 (0.58, 1 .0)' 








Halothane 


0.82 (0.26, 2.6)' 








Desflurane 


3.9 (2.1, 7.2)^ 









p-value from likelihood ratio test, ASA: American Society of Anesthesiologist; CI: Confidence interval; ETT-C: Endotracheal tube with controlled ventilation; LMA-A: 
Laryngeal mask airway with assisted ventilation; NDMR: Non-depolarizing muscle relaxant; URI: Upper respiratory tract infection, ^remote sites including radiology 
service and cardiac catheterization, hazard ratios within columns and variables that have no superscript {^^) in common differ significantly at p < 0.05 (Wald test). 



with IRE in our study. Previous studies did not distin- 
guish the effects in different periods [3-8]. The factors 
that appeared to have different effects in each period of 
anesthesia in our study were age, history of respiratory 
disease and ASA classification for any IRE, and intubating 
agent for laryngospasm. Patient position and anesthesia- 
related risk factors for IRE i.e. airway device with mode of 
ventilation, using succinylcholine for intubation and using 
desflurane anesthesia were somewhat surprising findings 



best described by analysing each IRE separately as was 
done in this study. 

Desaturation was the most common IRE (54%). A 
high incidence density of desaturation (44.9 per 100,000 
person-minutes) was seen during the induction period 
of anesthesia, a result supported by Bunchungmongkol et al. 
[12] and Hallan et al. [13]. Light anesthesia with ETT placed 
during the maintenance period (31%) was the most com- 
mon cause of desaturation, which might be precipitated 
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Table 3 Multivariate Cox regression of desaturation based on period of anesthesia (n = 42,055) 



Variable All periods (235) Induction (90) Maintenance (124) Emergence (21) 

Adjusted hazard ratio Adjusted hazard ratio Adjusted hazard ratio Adjusted hazard ratio 
(95% CI) p-value (95% CI) p-value (95% CI) p-value (95% CI) p-value 



Age (ref= > 6) 
1-6 


I'' < 0.001 
2.3 (1.6, 3.3)^ 


- - - 


< 1 


3.7 (2.5, 5.6)' 




Respiratory disease (ref = no) 




l'< 0.001 1'0.16 1'0.6 


URI 


- 


1 .1 (0.45, 2.9)' 1 .7 (0.89, 3.2)' 1 .5 (0.34, 6.8)' 


Pulmonary disease 




3.9 (2.4, 6.3)^ 1 .6 (0.84, 2.9)' 1 .9 (0.54, 6.5)' 


Site of operation (ref = superficial) 


l'< 0.001 




Broncho/laryngoscopy 


4.1 (2.1, 8.3)^ 




Ear-nose-throat 


1 .7 (0.87, 3.4)' 




Eye 


1.0 (0.49, 2.1)' 


- - - 


Thoracic 


1 .3 (0.62, 2.8)' 




Abdomen 


1 .3 (0.70, 2.4)' 




Spine & extremity 


1 .2 (0.63, 2.4)' 




Remote* 


1 .7 (0.64, 4.3)' 




Position (ref = supine) 


1' 0.01 




Lateral 


2.0 (1.2, 3.3)^ 


- - - 


Lithotomy 


0.42 (0.15, 1.2)' 




Prone 


0.77 (0.28, 2.1)' 




ASA classification (ref= 1) 




l'< 0.001 1'0.75 1'0.93 


2 


- 


1 .9 (0.92, 3.9)' 1 .1 (0.66,1 .7)' 1 .1 (0.56, 2.2)' 


3 




4.4 (2.0, 9.6)^ 1 .3 (0.70, 2.3)' 0.90 (0.25, 3.1 )' 


Airway device (ref = E^-C) 


1' 0.005 




Jet ventilation-C 


2.2 (0.89, 5.3)' 




Facemask-A 


2.5 (1.3, 4.7)^ 




LMA-A 


2.0 (1.1, 3.6)^ 




Intubating agent (ref=NDMR) 


l'< 0.001 




Succinylcholine 


2.0 (1.4, 2.8)^ 




Sevoflurane 


1 .2 (0.62, 2.4)' 




Inhalation agent (ref = sevoflurane) 


l'< 0.001 




Isoflurane 


0.79 (0.57, 1.1)' 




Halothane 


1 .2 (0.40, 3.8)' 




Desflurane 


5.7 (3.1, 10.7)^ 




p-value from likelihood ratio test, ASA: American Society of Anesthesiologist; CI: Confidence Interval; ETT-C: Endotracheal tube with controlled ventilation; LMA-A: 
Laryngeal mask airway with assisted ventilation; NDMR: Non-depolarlzing muscle relaxant; URI: Upper respiratory tract Infection, ^remote sites Including radiology 
service and cardiac catheterization, hazard ratios within columns and variables that have no superscript (^'^'^) In common differ significantly at p < 0.05 (Wald test). 


by positioning or surgical stimulation during the first precipitate laryngospasm at the time of surgical stimulation. 
15 minutes. Being in the lateral position increased the Assisted ventilation via facemask or LMA increased the 
hazard of desaturation two-fold compared to the supine hazard rate of laryngospasm more than 12-fold compared 
position. Laterality could affect ventilation and perfusion with controlled ventilation via ETT. However, our result 
matching and easily lead to desaturation aggravated by contrasts with the study of Flick and Bordet et al. [14,15], 
light anesthesia. in which the incidence of laryngospasm was found to be 
The second most common cause of desaturation was only slightly higher in children whose airway was managed 
laryngospasm (26%). A high rate of laryngospasm was seen with LMA (1.7%) compared to those with ETT (1.1%) and 
in the first 30 minutes when airway stimulation might actually lower in some other studies [8,16,17]. The true 
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Table 4 Multivariate Cox regression of wheezing based on period of anesthesia (n = 42,055) 



Variable All periods (101) Induction (48) Maintenance (41) Emergence (12) 

Adjusted hazard ratio Adjusted hazard ratio Adjusted hazard ratio Adjusted hazard ratio 
(95% CI) p-value (95% CI) p-value (95% CI) p-value (95% CI) p-value 



Respiratory disease (ref = no) 


1' < 0.001 


URI 


3.5 (2.0, 62t 


Pulmonary disease 


2.9 (1.7, 5.1)^ 


ASA classification (ref= 1) 


1^0.014 


2 


2.0 (0.99, 3.9)' 


3 


2.9 (1.4, 6.1)^ 



p-value from likelihood ratio test, ASA: American Society of Anesthesiologist; CI: Confidence interval; URI: Upper respiratory tract infection, hazard ratios within 
columns and variables that have no superscript {^^) in common differ significantly at p < 0.05 (Wald test). 



explanation why our result was different compared with 
other studies is unknown. However, a high rate of laryn- 
gospasm at the time of surgical stimulation could possibly 
occur in inadequately anesthetized children whose airway 
was managed by spontaneous breathing with assisted 
ventilation via facemask or LMA rather than those receiv- 
ing NMBA and controlled ventilation via ETT, Moreover, 
desaturation is a usual consequence of laryngospasm and 



could therefore be misdiagnosed as desaturation alone, 
especially when an ETT is in place. 

Use of desflurane incurred a hazard rate of desaturation 
and laryngospasm 5- and 10-fold higher, respectively, com- 
pared to the use of other inhalation agents, a result sup- 
ported by a study by Lerman et al. [18]. The mechanism of 
desflurane inducing airway contraction by activating transi- 
ent receptor potential-Al during spontaneous breathing 



Table 5 Multivariate Cox regression of laryngospasm based on period of anesthesia (n = 42,055) 


Variable 


All periods (75) 


Induction (27) 


Maintenance (38) 


Emergence (10) 




Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Adjusted hazard ratio 
(95% CI) p-value 


Age (ref= > 6) 




1' 0.08 


1' 0.004 


1' 0.002 


1-6 




1.4 (0.57, 3.6)' 


3.4 (1.5, 7.8)^ 


8.5 (1.8, 40.3)^ 


< 1 




3.7 (1.2, 11.3)^ 


3.4 (0.99, 1 1 .6)' 


27.4 (5.2, 143.2)' 


Respiratory disease (ref = no) 


1' 0.02 




s 




URI 


1 .7 (0.78, 3.7)' 








Pulmonary disease 


3.1 (1.5,6.4)^ 








ASA classification (ref= 1) 


1' 0.23 








2 


0.82 (0.47, 1 .4)' 








3 


0.39 (0.12, 1.3)' 








Airway device (ref = E^-C) 


l'< 0.001 








Jet ventilation-C 










Facemask-A 


18.1 (6.4, 51.4)^ 








LMA-A 


12.5 (4.6, 33.9)^ 








Intubating agent (ref = NDMR) 




l'< 0.001 


1' 0.67 




Succinylcholine 




13.5 (3.9, 46.3)^ 


1 .4 (0.28, 6.9)' 




Sevoflurane 




1.9 (0.18, 19.5)' 


2.8 (0.6, 13.0)' 




Inhalation agent (ref = sevoflurane) 


1' < 0.001 








Isoflurane 


0.75 (0.38, 1 .5)' 








Halothane 


1 .6 (0.37, 6.5)' 








Desflurane 


11.0 (5.1, 23.9)^ 









p-value from likelihood ratio test, ASA: American Society of Anesthesiologist; CI: Confidence interval; ETT-C: Endotracheal tube with controlled ventilation; LMA-A: 
Laryngeal mask airway with assisted ventilation; NDMR: Non-depolarizing muscle relaxant; URI: Upper respiratory tract infection, hazard ratios within columns and 
variables that have no superscript {^^^) in common differ significantly at p < 0.05. 
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Table 6 Multivariate Cox regression of upper airway obstruction in combination with reintubation based on period of 
anesthesia (n = 42,055) 



Variable 



All periods (23) 

Adjusted hazard ratio 
(95% CI) p-value 



Induction (1) 

Adjusted hazard ratio 
(95% CI) p-value 



Maintenance (15) 

Adjusted hazard ratio 
(95% CI) p-value 



Emergence (7) 

Adjusted hazard ratio 
(95% CI) p-value 



Age (ref= > 6) 

1-6 

< 1 

Site of operation (ref= superficial) 1^ < 0.001 

Bronclio/laryngoscopy 22.8 (2.6, 201 .8)^ 

Ear-nose-throat 4.7 (0.57, 38.3)' 
Eye 

Thoracic 0.83 (0.05, 13.3)' 

Abdomen 1.2 (0.13, 9.9)' 

Spine & extremity 1.4 (0.14, 13.7)' 
Remote* 



1' 0.01 

3.7 (0.90, 14.6)' 
9.2 (2.2, 39.0)^ 



1' 0.005 

4.7 (0.90, 23.3)' 

18.5 (3.7, 93.6)^ 



p-value from likelihood ratio test, CI: Confidence interval, ^remote sites including radiology service and cardiac catheterization, hazard ratios within columns and 
variables that have no superscript {^^) in common differ significantly at p < 0.05 (Wald test). 



with assisted ventilation via facemask or LMA makes lar- 
yngospasm more likely than when using airway protection 
with ETT [19]. 

Succinylcholine is used for ETT intubation in our 
institute in high risk children such as rapid sequence 
induction or difficult airway. ASA classification 3 and 
having respiratory disease were often considered as indi- 
cating high risk and were put into the Cox regression ana- 
lysis of desaturation and laryngospasm even though they 
were not significant predictors. The use of succinylcholine 
for intubation increased the hazard rate of desaturation 
two-fold higher compared to use of a non-depolarizing 
muscle relaxant, whereas the use of succinylcholine at 
induction period (first 15 minutes) incurred a hazard rate 
of laryngospasm much higher (13.5-fold) after adjusting 
for ASA classification and having respiratory disease. 
However, the temporal relationship between succinylcho- 
line and laryngospasm could not be determined because 
succinylcholine can be used as a treatment for laryngos- 
pasm [20] or for children requiring rapid sequence in- 
tubation [21]. Thus, the direction of causation should 
be interpreted with caution. 

In three multivariate Cox regression analyses, i.e. 
models having IRE, desaturation and laryngospasm 
as the outcome, the influence of airway device was 
consistent. ETT protected all three adverse events. 
The parameter of mechanical ventilation set by stand- 
ard operating procedures could not be estimated sep- 
arately because the mode of ventilation was highly 
linked to airway device i.e. children with ETT using 
NMBA received controlled ventilation whereas chil- 
dren who breathed spontaneously through facemask 
or LMA received assisted ventilation. It was not pos- 
sible to determine whether aiway device or mode of 



ventilation had a stronger effect on adverse events 
from our data. 

Wheezing or bronchospasm (22%) was the third most 
common cause of desaturation in our study. A high 
incidence density of wheezing or bronchospasm 
(23.9 per 100,000 person-minutes) was seen during the in- 
duction period of anesthesia. High risk children including 
having URI or pulmonary disease and ASA classification 3 
were the risk factors for wheezing in all periods in our 
study regardless of age and airway management, a finding 
contrary to other studies [8,16]. 

The incidence density of UAO in combination with 
reintubation was highest during the emergence period of 
anesthesia (3.4 per 100,000 person-minutes). UAO in 
combination with reintubation occurred mostly during 
60-180 minutes of anesthesia when the operation was 
finished and the airway device was taken out. The main 
risk factors were age < 1 years and airway surgery, a result 
supported by Ing et al. [22] and Rujirojindakul et al. [23] 
Anesthesia-related risk factors could not be determined 
because of the low number (23) of events. 

Strengths and limitations 

The strength of this study was its use of time-to-event 
analysis, which takes into account the changing risk of 
IRE over the duration of GA and incidence density per 
person-time, which has not been done in previous stud- 
ies in children. Because of using routinely recorded data, 
a limitation was that some patient-related risk factors, 
such as history of parental smoking or passive smoking 
and obstructive sleep apnea syndrome, were not available 
in the assessment. Although the historical cohort design 
may result in information bias from under-reporting, 
which can underestimate the true incidence [24], we 
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managed to determine statistically significant surgery- 
and anesthesia-related risk factors including those with 
time-varying effects. Moreover, the accuracy of the 
study in terms of internal validity was quite high due 
to the absence of false positives, the low rate of missing 
data (1.4%) and the large sample size. Because of the 
high internal validity, the generalizability to other child 
populations should be persuasive even though a single 
hospital was used to recruit subjects [25]. 

Conclusions 

The risk of most forms of IRE except UAO and of any 
IRE was highest during the induction and early mainten- 
ance periods of anesthesia (first 30 minutes of anesthesia). 
The common risk factors for IRE at any time were airway 
surgery, desflurane anesthesia and airway device with mode 
of ventilation. The effect- varying risk factors of any IRE 
were age, history of pulmonary disease and ASA classifica- 
tion. Therefore, anesthesiologists should pay more attention 
during the induction and early maintenance periods 
especially when certain airway devices incorporating 
assisted ventilation or desflurane are used. 
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Additional file 1: Univariate Kaplan Meier analysis of intraoperative 
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